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The feasibility of complete photoionization experiments, in which the full set of photoionization 
matrix elements are determined, using multiphoton ionization schemes with polarization-shaped 
pulses has recently been demonstrated [Hockett et. ah, Phys. Rev. Lett. 112, 223001 (2014)]. Here 
we extend on our previous work to discuss further details of the numerics and analysis methodology 
utilised, and compare the results directly to new tomographic photoelectron measurements, which 
provide a more sensitive test of the validity of the results. In so doing we discuss in detail the 
physics of the photoionziation process, and suggest various avenues and prospects for this coherent 
multiplexing methodology. 


I. INTRODUCTION 

The aim of “complete” photoionization studies is the de¬ 
termination of the amplitudes and phases of the ionization 
matrix elements, which constitute a fundamental descrip¬ 
tion of an ionization event Bill. The matrix elements 
define the coupling of the initial state to the final com¬ 
pound state, comprised of an ion and free electron. In 
the dipole limit, this matrix element can be very generally 
defined as 4 ^+|/x.E| 4^). Here 4^ is the initial wavefunc- 
tion of the system, the photoion, the photoelectron, 
fi the dipole operator and E the electric field. By express¬ 
ing the continuum wavefunction as a set of partial-waves^ 
corresponding to different continuum angular momentum 
states, the ionization matrix element can be decomposed 
into various geometric and radial components, and the set 
of amplitudes and phases of these components constitutes 
a complete description of the ionization event. In order to 
determine these matrix elements from experimental data, 
an observable sensitive to the relative phases of the partial- 
waves is required, and such an interferometric observable is 
found in the photoelectron angular distributions (PADs), 
which are angular interference patterns dependent on the 
composition of 

A range of experiments have been performed in order 
to provide such complete descriptions of photoionization 
for a number of atomic and molecular systems. The key 
concern in such experiments is the level of detail required 
in order to undertake the relatively complex analysis pro¬ 
cedure. Typically the angular (or geometric) part of the 
matrix elements can be calculated analytically [3], leaving 
only the energy-dependent radial (or dynamical) compon- 
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ents to be determined from the experimental data. The 
determination of these components involves fitting experi¬ 
mental data with the specific ionization formalism for the 
ionization event under study. Since, in general, there may 
be many partial-waves and the composition of is not 
usually known a priori, a large experimental dataset is re¬ 
quired for this procedure. In order to obtain a sufficient 
dataset, experimental data is obtained for a range of geo¬ 
metric parameters, for example by varying the polarization 
state and polarization geometry mg or, for molecules, the 
rotational state or axis distribution mig, or via molecu¬ 
lar frame measurements [HHis]. Since the dynamical para¬ 
meters are invariant to these geometric changes, a dataset 
of sufficient information content to determine these para¬ 
meters may be obtained in this way. 

Recently, we demonstrated a new type of measure¬ 
ment and analysis methodology for complete experiments 
m- This method can be considered as time-domain 
polarization-multiplexing. In this case, a multiphoton ion¬ 
ization scheme with a moderately intense, ultrafast laser 
pulse was employed to ionize potassium atoms. The result¬ 
ing light-matter interaction can be understood as an intra¬ 
pulse two-step processes, in which electronic population 
transfer is driven by the laser field (i.e. Rabi oscillations), 
and the excited state population created can subsequently 
be ionized via 2-photon absorption. In this case, the popu¬ 
lation dynamics and the ionization dynamics are dependent 
on the properties of the laser pulse, as well as the physical 
properties of the system which ultimately determine the 
matrix elements. In this scheme, changing the polarization 
of the pulse corresponds to changing the geometric para¬ 
meters of the ionization, as described above. In the simplest 
case a single or pure polarization state is employed, and the 
geometric parameters are time-invariant. More generally, 
via the use of a polarization-shaped pulse, the geometric 
parameters can be changed in a time-dependent manner. 
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Since the dynamics and ionization all occur within a single 
laser pulse, the process is fully-coherent, and the final, time- 
integrated, photoelectron measurement can be considered 
as a time-domain multiplexed measurement of the set of 
(instantaneous) polarization states explored by the shaped- 
pulse. 

Here we discuss further details of the work presented in 
ref. [17], with a focus on extending the details of the the¬ 
ory presented therein, in particular the numerical details of 
the fitting procedure and a discussion of the benefits and 
limitations of this approach. We further present detailed 
comparison of our results with new maximum information 
photoelectron measurements, utilizing a tomographic pro¬ 
cedure for the measurement of 3D photoelectron distribu¬ 
tions and detailed analysis, allowing for a quantitative com¬ 
parison of the predicted PADs and experimental PADs as 
a function of polarization geometry (further details of the 
maximum information measurements can be found in ref. 
m)- Finally, the possibilities of extending this treatment 
to different classes of ionization is explored, with a partic¬ 
ular emphasis on molecular ionization problems. 

II. INTRAPULSE DYNAMICS MULTIPHOTON 

IONIZATION WITH POLARIZATION-SHAPED 
PULSES 

Here we detail the various steps involved in the treatment 
of the 3-photon scheme detailed above. For completeness 
we include all aspects of our treatment. 
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A. Electric field 

The electric field as a function of time is described as: 


E{t) = (1) 


Figure 1. Time-dependent dynamics for an elliptically polarized 
laser pulse, (a) Laser field envelope for an elliptically polarized 
pulse, defined by = 0.5 rad (eqn. [^. (b) Bound state pop¬ 
ulations (eqn. (c) Instantaneous continuum populations 

dif^rnf(k,t) (eqn. TTl) and (d) cumulative continuum population 
(eqn. [l^ . 


where Eq is the field strength, the pulse envelope is Gaus¬ 
sian with temporal width parameter r and uj is the carrier 
(angular) frequency. Using the notation of ref. [19] the 
spectral content of the pulse is given by: 

E{n) = E{Eit)} ( 2 ) 



The field can also be expressed in terms of a spherical 
basis, i.e. left and right circularly polarized components: 


where T represents a Fourier Transform. 

Polarization shaped pulses are described as in ref. m, 
by assuming initially identical and in-phase (x, y) field com¬ 
ponents, then applying a spectral phase shift. Hence, a 
field described by two Cartesian components with inde¬ 
pendent spectral phases (but identical spectral content) can 
be defined as: 


f E,{n) 

V Eyin) 


= E{Q) 


^i(j)y{Q) J 


resulting in the time-domain components: 


(3) 


-£^l(^) a _ ^ f Ex{t) iEy(t) \ / X 

Enit) J - ^[Ex{t)EiEy{t) J 

This final form was used in the calculations herein, since it 
physically describes the instantaneous pulse angular mo¬ 
mentum, in terms of the projection of the photon mo¬ 
mentum onto the propagation axis, where L equates to 
m = +1 and R to m = — 1 states. This form can therefore 
be directly interpreted in terms of the allowed Am of both 
bound-bound and bound-free transitions - this is discussed 
further below. Note that this form implies that the light 
propagates along the 2 :-axis, and the lab. frame angular 
momentum m is defined relative to this propagation axis. 
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Figure 2. Time-dependent dynamics for a polarization-shaped 
laser pulse, (a) Laser field envelope for an elliptically polarized 
pulse, defined by (l)y = —tt rad for the blue half of the pulse (eqn. 
[^. (b) Bound state populations (eqn. [^. (c) Instantaneous 
continuum populations t) (eqn. rTTl) and (d) cumulative 

continuum population (eqn. |12[). 


B. Non-perturbative laser-atom interaction 


The strong laser field drives Rabi oscillations in the atom, 
coupling electronic states |n, /, m). In the case of potassium 
atoms, as detailed in ref. [20], the initial population is 
in the 45 state and the laser frequency is near resonant 
with the |4, 0, 0) ^ |4,1, m) transition, hence single photon 
absorption populates the 4p manifold, while a strong laser 
field will drive Rabi cycling between the 45 and 4p states. 
The allowed values of m depend on the polarization state 
of the light. 

The population dynamics during the laser pulse, de¬ 
scribed in the spherical basis of eqn. are given by the 
time-dependent Schrodinger equation: 


d ( \ ( ^ \ ( S{t) 

-f. P+i{t) =M <5+1 0 

0 <5_i J ^ 

(6) 

where 5(t), p+i(t) and P-i{t) are the state vector com¬ 
ponents for the |4,0,0) and |4,1,±1) states, ^L/nit) = 
where are the transition amplitudes, 

and 6±i represent the detuning of the laser from the res¬ 
onant frequency of the transition. Here it is clear that the 
L and R components of the electric field drive transitions 
with Am = +1 and Am = — 1 respectively; this is simply 
the consequence of the conservation of angular momentum 
since the light carries I = 1 unit of angular momentum, 
with lab. frame projection m = 1 for El and m = — 1 for 
Er. In this sense the (instantaneous) helicity of the electric 
field is directly imprinted on the atomic ensemble. 

Here h and Eq are both set to unity for simplicity; pl/r is 
also set to unity, i.e. equal probability of transitions to both 
|4, l,m) states, and 6±i = 0.05 rad/fs. For determination 
of PADs these simplifications are acceptable as only the 
relative population of m = ±1 states will affect the angular 
distribution, and these populations are dependent only on 
the driving laser field polarization. 

C. Perturbative two-photon ionization PADs 

In the perturbative regime, the dipole transition amp¬ 
litude for a transition from a bound state \ni^li^mi) to a 
continuum state |k; If^mf) is given by the dipole matrix 
elements: 


di^f{kR) = (k; lfmf\fi^f.E{t)\nilimi) (7) 

« RZif{k)Eg{t){lfmf, Iqlkmi) (8) 

where is the dipole operator; is the radial part 

of the matrix element, which is dependent on the mag¬ 
nitude of the photoelectron wavevector k, the principal 
quantum number of the initial state n and the electronic or¬ 
bital angular momentum /, but assumed to be independent 
of mi and mf; {lfmf,lq\limi) is a Clebsch-Gordan coeffi¬ 
cient which describes the angular momentum coupling for 
single photon absorption, with q = ±1 for the L and R 
components of the laser field respectively. This treatment 
corresponds to a single active electron picture, in which the 
final state is a pure continuum state, i.e. the photoion is 
neglected and there is no angular momentum transfer to 
core. Spin is also neglected. This treatment is sufficient 
for the potassium atom case discussed herein; extension to 
more complex coupling schemes is discussed in sect. [V| 
Under these assumptions, the angular part of both 
bound-bound and bound-free transitions are described by 
matrix elements of the same form. Using these dipole mat¬ 
rix elements, two-photon ionization to a single final state 
from an initial state via a virtual one- 

photon state can then be written as: 
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difTTif — di^y(^k^t^dy^f(^k^t^ 


— ^ ^ -^Z^Z j (^) (^/^/5 l^v'^v) {^v'^v •> ^Q\^i'^i)Xni,li,mi{t^ 

Tli , li , TTLi 

Tly , ly , TTly 

q^q' 


(9) 

( 10 ) 


This form shows the general case, with summation over 
all initial states \niji,mi) weighted by their populations 
Xni,ii,mi{t)' Although the bound-free matrix element is la¬ 
belled with quantum number in practice this is un¬ 
assigned and will correspond to a quasi-continuum of vir¬ 
tual states within the laser bandwidth, so is dropped in 


the following. El In this treatment all energy dependence 
is contained in the R{k) radial integrals. For the potassium 
case considered here, a slightly simplified form can be writ¬ 
ten since only the 4p levels contribute to the ionization, 
hence rii = 4, and the time-dependent populations are given 
by Prrii (t) as defined in eqn. 


E f ik)Eq' (j:,^ (^l fTTif ^ \q \lyTny)Ri^^j^^(^k^Eq(t^(^ly7Tiy^ lq\li7Tii)prfji.(t^ (H) 

1 

ly , TTly 

q^q' 

Integrating over t yields: 

^Z/m/(^) di^y(^kR^dy^f(^kR^dt ^ ^ \lv'^v)Ej^^j^^(^k^(^lyTfiy^\q\liTfii^Eqf(t^Eq(t^pyiri^(J^^ 

li , ITli , ly , n^v 

q,q' 

( 12 ) 

The observed photoelectron yield as a function of angle, the PAD, for a small energy range dk over which we assume 
the R{k) are constant, is then given by the coherent square over all final (photoelectron) states: 


I{e,4>-, k) = fdk E {0, 4>)d;,^, {k)Y^^^, {0, 4>) 

lyrrif 

lym'f 


(13) 


This treatment is very similar to that given in ref. 
m, with the Clebsch-Gordan coefficients equivalent to the 
<^z,m;Z',m' parameters and the di^mf similar to the cz,m- 
The main difference is that all are accounted for, 

hence the explicit inclusion of the radial elements Rii{k). 
The radial matrix elements defined here are assumed to 
be complex, and include both the scattering phase 
and the geometric phase factor which usually appear in 
the definition of the photoelectron wavefunction [22]. The 
amplitudes and phases of these parameters constitute the 
unknowns which are sought in “complete” photoionization 
studies and, physically, define the scattering of the outgoing 


photoelectron from the nascent ion core. 

The PAD can also be described by a generic expansion 
in spherical harmonics with expansion coefficients 
termed anisotropy parameters, where: 


I{0, cP;k) = J2 l3L,M{k)YL,M{0, 4>) (14) 

L,M 

In general the PL,M{k) provide a compact way to express 
the PADs, and allowed values are constrained by symmetry 
[U [53] . This expansion can be considered as indicating the 
information content of a given distribution, and the res¬ 
ultant multipole moments L, M are related to the partial 
wave expansion of eqn. 14 by [9]: 


h,M 


/(2//-1-I)(2/^-1-I)(2L + I) / ij^ I'^ L 

^ V ^ 0 0 0 

/j, rrif 


If I'f L 

nif —m'f M 


difmf{k)d*,^^,^ 


(k) 


(15) 
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Further exploration of the information content of PADs 
for the case of tomographic 3D photoelectron measure¬ 
ments can be found in ref. m 


D. Pure and shaped laser pulse dynamics 


In order to illustrate the theory detailed above, figures 
and give the details of two example calculations, for an 
elliptically polarized pulse and a fully polarization-shaped 
pulse respectively. In both cases the panels illustrate, from 
top to bottom, the envelope of the laser field and L, R 
components, as defined by eqn. the population dynam¬ 
ics driven by the laser field, in terms of the state vector 
components s(t), p+i(t) and P-i{t) defined in eqn. the 
instantaneous continuum populations, as defined by eqn. 
TT] and making use of the previously determined p hot oion- 
ization matrix elements Ru (see ref. m and sect. pTO; the 
cumulative continuum populations, as defined by eqn. El 

Both examples provide insight into the dynamics of the 
ionization process, and it is clear how the L and R compon¬ 
ents of the laser field drive both the bound-state population 
dynamics, and the instantaneous continuum contributions. 
Since, in this model, the two steps are decoupled, and the 
ionization is assumed to be perturbative, there is no deple¬ 
tion in the bound state populations. The ionization step 
does, however, follow the bound state dynamics since the 
instantaneous population defines which continuum states 
can be accessed, and their relative weighting. Thus, the 
instantaneous continuum dynamics follow the bound-state 
dynamics. Furthermore, since there are no continuum elec¬ 
tron dynamics in this model (i.e. no laser-continuum coup¬ 
ling, or electron-ion recombination), the final continuum is 
simply the sum over the instantaneous continuum contri¬ 
butions (eqn. and builds up coherently over the pulse 
envelope. The resultant PAD, eqn. thus depends both 
on the final continuum populations, as well as the accumu¬ 
lated phase for each state. 

In the case of a “pure” polarization state (fig. 0 , in 
this example an elliptically polarized light field defined by 
4>y = 0.5 rad., there is essentially no dynamic contribution 
to the final result since the relative continuum contribution 
is time-independent. In the language used previously, the 
geometric contribution to the ionization is time-invariant. 
However, in the case of a polarization-shaped pulse (fig. 
[^, where the relative L and R components do vary sig¬ 
nificantly over the pulse, the intra-pulse dynamics play a 
key role in defining the final continuum wavefunction. It 
is this dependence that makes the final PAD particularly 
sensitive to the pulse shape, as well as the ionization mat¬ 
rix elements. While the two cases are formally identical, 
there is clearly no polarization multiplexing in the pure 
case, since the polarization state is time-invariant. In the 
polarization-shaped case, the information content is greatly 
increased since the final result arises from coherent addition 
over all instantaneous polarization states, thus contains ad¬ 
ditional information relative to a pure case. (Further ex¬ 
amples of polarization-shaped pulses and resultant PADs 
can be found in ref. ini-) 


III. PHOTOELECTRON IMAGE GENERATION 
FITTING 

In this section we outline salient details of the numer¬ 
ics used in applying the above theory to the generation of 
photoelectron momentum distributions which can be com¬ 
pared with experimental data. In the context of complete 
photoionization experiments, the use of these momentum 
distributions to generate 2D photoelectron images and fit 
experimental data is described. 


A. Photoelectron momentum distributions 

The theory detailed above provides a definition of the 
photoelectron yield as a function of time, energy and angle, 
most compactly defined by the /3 l,m(^) parameters, but ul¬ 
timately depending on the underlying laser and target prop¬ 
erties. The generation of theoretical, time-integrated, pho¬ 
toelectron momentum distributions from these parameters 
simply involves the population of a 3D grid ( 6 >, 0, k) with 
the relevant basis set expansion in spherical harmonics as a 
function of energy, as defined in eqn. (the radial aspect 
of this expansion is discussed below). 

The volumetric data defined in this way is equivalent to 
the experimental data recorded in a 3D imaging experi¬ 
ment, examples of such experiments are direct 3D imaging 
via techniques with high temporal and spatial resolution 
(for instance refs. [24]-[26] and references therein), or indir¬ 
ect methods based on tomography in which 3D distribu¬ 
tions are reconstructed from a set of 2D projections [2M9] 
(see also ref. [H]). For comparison with 2D imaging data, 
further integration along a spatial dimension is additionally 
required in order to project the volumetric data onto a 2D 
plane. We note that in both imaging experiments and the 
numerics applied here, this summation is treated incoher¬ 
ently. Physically, this corresponds to a loss of photoelectron 
coherence before or at the detector, effectively long after 
the coherent quantum mechanical scattering event which 
determines the momentum distribution (PADs and energy 
spectrum) [30l|3T]. Since the range of the initial scattering 
event is microscopic, while photoelectron propagation and 
detection is macroscopic and often involves the application 
of external fields and, ultimately, discrete particle counting, 
this is a physically reasonable assumption. 

In the results shown in paper I we additionally as¬ 
sumed that the radial dependence of the ionization mat¬ 
rix elements over the span of the main spectral feature 
(~200 meV) was negligible, and that the details of the ra¬ 
dial distribution could be simplified to a Gaussian energy 
spread with no phase contribution. This allowed for the 
momentum data generation and fitting to be simplified, 
and the radial distribution given by a Gaussian (defined in 
energy-space): 


G(k) = (16) 

V27r7 

where /q is the intensity, 7 the width, E{k) and E{ko) define 
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Figure 3. Illustration of computed PADs and 2D photoelectron images as a function of pulse polarization defined by (j)y^ propagating 
in the ^-direction. The top row shows PADs in polar form, /(^, 0; /c), the bottom row shows the same distributions projected onto 
a velocity iso-sphere (single /c), and 2D image plane projections I2d{02d, k) assuming a Gaussian energy spectrum. 


the radial coordinate and the peak centre in energy-space. 
The final 3D momentum distribution is then defined by: 

I{0,^,k)=G{k)J2l^iMYL,M{O,ct>) (17) 

L,M 

Where the /3^ ^ include the superscript to denote that 
these parameters are generally dependent on k, as in eqn. 

but are here taken to be constant over the range of 
k spanned by the Gaussian envelope G{k). Finally, it is of 
note that more generally the Gaussian assumed here should 
be replaced by an accurate energy spectrum, this point is 
discussed in sect. El 

The 2D images obtained by integration of the volumetric 
distribution function are then given as: 

hD{02D,k) = J l{0,(l),k)du (18) 

u 

where u defines the domain of integration (with integration 
over the Cartesian X^Y oi Z directions for the correspond¬ 
ing {y,z), {x^z) or (x^y) image planes respectively), and 
O 2 D is defined in the image plane. 

Figure illustrates the computed PADs, obtained us¬ 
ing the matrix elements of ref. [H], for four polarization 
states of the electric field. The top row shows the PADs 
in spherical polar form, as defined by eqn. while the 
bottom row shows the same PADs projected onto spherical 
surfaces. This is how the distributions appear in velocity 
space, as the angle-dependent photoelectron flux for each 
k. The 2D projections show the same angular distribu¬ 
tions, combined with a Gaussian energy spectrum as per 
eqn. Eli and projected onto 2D Cartesian planes. These 


image planes simulate velocity map imaging data, and il¬ 
lustrate how the experimental results will depend on both 
the native details of the PAD and the details of the pro¬ 
jection geometry. In this case, the laser propagates along 
the Z-axis, and the polarization is defined in the {X^Y) 
plane, so experimental images will correspond to the image 
planes (x, z) or (^, z) (since images cannot be obtained in 
the propagation direction in a standard VMI experiment), 
and the precise details will further depend on the rotation 
of the distribution about the Z-axis. (Further details of 2D 
and 3D imaging, geometry considerations and information 
content, can be found in ref. EH-) 


B. Fitting methodology 

As discussed above, within the framework developed 
herein the radial matrix elements are the only unknown 
quantities. With a sufficient experimental dataset one can 
therefore hope to obtain these matrix elements via a fit 
to the data. In this case the results of such a fit have 
already been presented in ref. m, and validated via good 
agreement with both the original 2D imaging data, and 
additional 3D data obtained via tomographic imaging ex¬ 
periments. We discuss here further details of the fitting 
methodology applied, since in general it is necessary to 
approach this complicated problem carefully. In partic¬ 
ular, we applied statistical analysis methodologies which 
were previously developed for energy-domain photoioniza¬ 
tion experiments [I1132]. 

In our procedure, the data from 2D measurements was 
compared with the calculated 2D images, as illustrated in 
figure and obtained as detailed above. The calculated 
images were then optimized via a fitting routine, with the 


7 




Figure 4. Example of fit statistics for Rp^d, (a) magnitude, (b) 
phase. 


radial matrix elements and image generation parameters as 
the free parameters for fitting. The criteria for the best fit 
was simply the minimization of the sum of least squares: 

{hD{02D,k) - lExpt.{^2D,k))‘^ (19) 

02D,k 

where I 2 D is the calculated distribution defined in eqn. 
and lExpt. the 2D experimental data. This methodology is 
completely general, and only relies on the underlying theor¬ 
etical framework correctly describing the physics inherent 
to the problem. However, the size of the hyperspace may 
be very large since it has dimensions equal to the number 
of free fitting parameters. The practical outcome of this 
is that the possibility of local minima in the hyperspace is 
significant, and the parameters obtained via such a proced¬ 
ure must be carefully evaluated and tested to confirm their 
veracity and robustness. 

In this particular case, the full calculation required 12 
parameters, consisting of the amplitudes and phases of the 
5 radial matrix elements Ru and 2 image generation para¬ 
meters (Gaussian centre and FWHM) [33]. Since absolute 




A(a.) 


Figure 5. ID cuts through the hypersurface for (a) mag¬ 
nitudes and (b) phases of the Ri^f matrix elements. Cuts 
are made by varying each parameter by Aaj, here given as 
a fractional variation from the best fit value (i.e. test value 
a'j = aj -h {aj x Aaj), and a' = aj for Aaj = 0), and the result¬ 
ing change in y^ is given as a percentage relative to the best fit 
value. 


phases cannot be determined, one phase is chosen to be 
a reference and set to zero, leaving 11 free fit parameters. 
Furthermore, the image generation parameters do not have 
a large influence on the final results, which are primarily 
sensitive to the angular coordinate, and could therefore be 
bounded quite tightly after some initial by eye optimiza¬ 
tion, thereby reducing the search-space of physical relev¬ 
ance to, effectively, 9 dimensions. In the fitting proced¬ 
ure the Ru were expressed in magnitude and phase form, 
Ru = where 0 < Ru < 1, -n < Su < tt, and 

Re^i^ = Ri^^i_^. Fitting was implemented with a stand¬ 
ard fitting algorithm, Matlab’s lsqcurvefit, based on a 
Trust-Region-Reflective least-squares method. The data¬ 
set for fitting consisted of four experimental images, each 
corresponding to a different pure polarization state of the 
light, similar to the states shown in figureThe laser pulse 
was modelled with r = 30 fs and four polarization states 
given by = 0 (linear polarization, ellipticity 5 = 0), 
4>y = TTjS^ 7T14: (elliptical polarization states, with ellipt- 
icities 5 ^0.2, 0.4) and (py = 7 r /2 (circular polarization, 
5 = 1). The ellipticities given here are defined as the ra¬ 
tio of the minor to major axes of the polarization ellipse, 
hence 5 = 0 for linearly polarized light and 5 = 1 for pure 




































































Figure 6. Calculated 3D distributions 1(0, cj), k) (first column) and sensitivity to changes in the phases and for (a) linear 
polarization and (b) elliptical polarization 0^, = 0.5 rad. 


circularly polarized light. Because the elliptical polariza¬ 
tion states may be slightly different from those obtained 
experimentally (via the use of a quarter-wave plate, see ref. 
m and refs, therein for further experimental details) the 
subsequent fitting was weighted towards the linear and cir¬ 
cular polarization results by an additional factor of two in 

In order to carefully test for local minima the hyperspace 
was repeatedly sampled using a Monte-Carlo approach, in 
which the fitting was repeated N-times with the seed val¬ 
ues for the fitting parameters randomized on each iteration. 
Statistical analysis of the fitted parameters derived from 
such repeated fits can be employed to probe the behaviour 
of the fitting algorithm, and also to gain information on 
how well the experimental data defines each fitted para¬ 
meter. Although it is non-trivial to visualize the full 
hypersurface, aspects can be probed by plotting histograms 
and correlation plots of the fitted parameters. A large scat¬ 
ter in the value of a given fit parameter over a range of fits 
to the same data suggests a poorly defined parameter; a 
consistent result meanwhile shows that a particular para¬ 
meter is well defined by the dataset. The experimental 
data can show different sensitivities to different parameters 
depending on the type of ionizing transitions present, be¬ 
cause different transitions will (according to the magnitude 
of the geometrical parameters and symmetry constraints) 
be more sensitive to certain partial-waves. Additionally, 
the presence of multiple-minima in the fit may be revealed 


by the presence of more than one feature in the histogram, 
reflecting more than one “best” fit result, while correlations 
appearing between supposedly uncorrelated parameters can 
indicate emergent behaviours in the high-dimensional space 
or - more prosaically - issues with the fitting methodology 
or coding. 

In this case we performed 300 fits, and the lowest was 
obtained on 4 of these fits, which we take to be the absolute 
minimum. The radial matrix elements ultimately found, as 
reported in ref. m, are given in table for reference, 
and discussed further below. Figure]^ gives an illustrative 
example of the fitting statistics, in this case showing cor¬ 
relation histograms between and the (a) magnitude and 
(b) phase of Rp^d- The plot shows the 40 fit results within 
5% of the lowest Interestingly, in this case many of the 
fits are bunched, with x^ ^ 300 (arb. units). This most 
likely reflects the presence of local minima as defined above, 
but may also be related to the convergence criteria set on 
the fitting algorithm which, in this case, was set to a lim¬ 
ited number of iterations in order to cap the computational 
time per fit and ensure a large seed-space for the search; 
in effect the large seed-space becomes part of the fitting 
criteria. Depending on the seed values, the overall conver¬ 
gence of the fit may be fast or slow, and the possibility of 
finding the global minima will vary depending on the start 
position in the 11-dimensional parameter-space, as well as 
the topography of this space and the details of the fitting 
algorithm. From the histograms of the bunched results. 
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d 
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0.14 (2) 

2(2) 

-2.08 (8) 
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f 

0.51 (9) 

26 (13) 

0.24 (7) 


Table L Fitted values for the relative transition matrix element 
magnitudes, \Rii\, and phases, Su. The square of the magnitudes 
is expressed as a percentage of the total transition amplitude, 
normalized to unity for each step. Uncertainties in the last digit 
are given in parentheses. * reference phase, set to zero during 
fitting. 


it is apparent that \Rp^d\ is somewhat well-defined at the 
larger with values mostly close to the best result, while 
the phase appears much less well-defined at this level. At 
lower the parameter-space is much sparser, with only 
a few parameter sets found, but they appear to converge 
on a single parameter set. These observations illustrate the 
difficulty in assessing best fit results without careful ana¬ 
lysis: in this case sampling of only a few fit results would 
potentially lead to a parameter set quite different from the 
global optimal found. Here statistical analysis, as well as 
further validation of the results against additional experi¬ 
mental data (see sect. IV), both serve to provide confidence 
that the absolute best fit, hence physically correct, results 
have been obtained. 


C. Robustness, uncertainties validation 

As well as statistically evaluating fit results, the beha¬ 
viour of x^ can also be more directly probed. In essence, 
this amounts to removing the black-box nature of the fit¬ 
ting algorithm by explicitly looking at the gradient and 
curvature of x^ as a function of the fitting parameters, 
rather than looking at only the final fitted results. Ad¬ 
ditionally, the curvature with respect to a given parameter 
can be used to provide uncertainty estimates on the fitted 
parameters [34] : 



( 20 ) 


where aj is the uncertainty in parameter aj. Equation 


sharper the response the better aj is defined by the data 
and hence the smaller the uncertainty. In practice this pro¬ 
cedure equates to varying each fitted parameter by Aaj, 
and evaluating x^ for this new parameter set, in order to 
map out ID cuts through the hypersurface. Uncertain¬ 
ties estimated in this manner were given in ref. and 

are provided again in table [T[ It is also of note that a similar, 
but not identical, procedure can be performed by refitting 
all other {n—1) parameters as a function of a test parameter 


20 relates the response of x^ lo a given parameter; the 


aj m- This procedure will also provide ID cuts through 
the hypersurface, but along the n-dimensional topography 
of the minimum. The drawback of this alternative pro¬ 
cedure is the necessity of performing many additional fits, 
which may be computationally expensive; for this reason it 
was not explored in this work. 

Figure shows ID cuts through the x^ hyperspace as 
defined above, for the magnitudes and phases of Ri^f. In 
this case it is clear that the sensitivity of x^ is good in most 
cases, with 10 % changes in aj (i.e. Aaj = ±0.1) typically 
leading to clear changes in x^; this is also reflected in the re¬ 
latively small uncertainties given in table hi In this case, 
a notable exception is |R 2 ^ 3 |, which is mucn less sensitive 
to Aaj for increases in magnitude. This is the magnitude of 
the /-wave channel, which dominates the ionization overall; 
consequently the final PAD is not very sensitive to small 
increases in the magnitude of this matrix element, although 
does remain very sensitive to decreases in magnitude, and 
its relative phase. In general x^ is somewhat less sensitive 
to the phases than the magnitudes, although the response 
is still significant. It is also of note that the ID cuts are 
not symmetric about Auj, reflecting the complicated to¬ 
pography of the x^ hypersurface, and the fact that it is 
dependent on the relative, rather than absolute, values of 
the matrix elements. 

A final, valuable test of the determined matrix elements 
is their predictive power, and the possibility of testing 
such predections against additional experimental results 
not used in the original extraction procedure. A consider¬ 
ation of the sensitivity of the determined matrix elements 
in these terms is a useful way of evaluating the results. In 
previous, energy domain, studies the (rotational) energy 
spectrum could be used to provide the additional, inde¬ 
pendent data against which the extracted matrix elements 
could be further verified m, and the possibility of us¬ 
ing different polarization geometries combined with tomo- 
graphically reconstructed PADs was also explored [29]. As 
noted above, and discussed briefly in ref. m , comparison 
of the current results with 3D photoelectron data obtained 
via tomography was also employed in this case. The com¬ 
parison with the experimental data is discussed in sect. 
IV B[ while the sensitivity of the computed 3D distribu¬ 
tions to changes in the matrix elements Aaj is discussed 
here. Figure provides some examples of this sensitivity 
for the variation of two different phases by 20 %, and for 
two different polarizations. Although the sensitivity of the 
3D distributions to these phases is inherent in the small 
uncertainties determined above, as well as the ability to 
successfully use a fitting methodology, it is nonetheless in¬ 
structive to visualize the sensitivity in this way. Here it 
is clear that, while both cases exhibit a sensitivity to the 
phase adjustments, the changes in the linearly polarized 
case are less significant. In this case, the width of the cent¬ 
ral bands increases slightly and, although this change still 
correlates with a change in the magnitude of 

this change means that it will only be revealed by careful 
quantitative analysis, and may not be obvious in a qualit¬ 
ative comparison. This conclusion becomes even stronger 
for 2D images of this distribution. In the elliptically polar- 



















10 


ized case the phase changes are manifested in an increase 
in flux and spread of the equatorial lobes of the distribu¬ 
tion, and are much more pronounced as compared to the 
linearly polarized case. While the sensitivity observed here 
merely confirms the earlier analysis, the investigation of 
the predicted distributions in this phenomenological man¬ 
ner provides additional insight into the fitting process, in 
particular the magnitude of changes which might be ex¬ 
pected in a given case and, hence, suggests possibilities for 
future experimental work, particularly in the more complex 
case of shaped pulses, as discussed in ref. m- 

Overall the methodology outlined here might be viewed 
as a pragmatic approach to complete experiments. Utiliz¬ 
ing a combination of fitting, Monte-Carlo sampling, direct 
exploration of the hyperspace and further validation of 
the results based on their predictive power, a careful valid¬ 
ation of their robustness and validity can be made for the 
case at hand. This is distinct from a more formal treatment, 
such as that discussed in Schmidtke et. al. m, wherein 
the fundamental limits of a fitting approach are derived. In 
the current work a comparison with the definitions given in 
that work have not been made, but the pragmatic meth¬ 
odology herein indicates that the extraction of the matrix 
elements is, in this case, reliable. An extension of this meth¬ 
odology, combined with a formal treatment, to investigate 
the additional possibilities in the polarization-multiplexed 
case remains for future work. 


IV. COMPARISONS WITH TOMOGRAPHIC 
DATA 

Here we focus on a detailed comparison of the results 
of ref. [T7] with additional experiments which provided 
full 3D data. These results, obtained using photoelectron 
tomography techniques (see ref. m for details), provide 
both a highly detailed volumetric data and a set of meas¬ 
urements at a different laser intensity. The former char¬ 
acteristic allows for a qualitative visual comparison of 3D 
distributions, which reveal details of the distributions which 
may be obscured in the 2D images, and the possibility of 
a full retrieval of the from the data, which is not 

possible for non-cylindrically symmetric 2D images and al¬ 
lows for a more quantitative comparison of experiment and 
theory. The use of different intensities 10^^ Wcm“^ for 
the tomographic data, as compared to 10^^ Wcm ^ for 
the 2D data) provides further evidence for the lack of any 
significant strong-held effects to the angular distributions 
in this case, and the veracity of our ionization model. 


A. Qualitative comparison 

A qualitative comparison of data for the three different 
polarization states measured is presented in hgure In 
this hgure the full 3D distributions are shown as nested 
iso-surface plots, and 2D images in the polarization plane 
are also shown. In this case, the experimental data has an 
additional high-energy feature in the radial distribution. 


arising from Antler-Townes splitting which becomes signi- 
hcant at higher intensities (see refs. uniEi, and ref. CHD- 
In this analysis only the main feature is of interest, and the 
tomographic distributions shown include a radial mask in 
order to remove the additional contributions and facilitate 
comparison over the main spectral feature. For the 2D im¬ 
ages (bottom row of hg. no radial mask is employed and, 
consequently, the experimental results show a broadening 
of the spectrum in the 2D images. Full details of the exper¬ 
imental data and tomographic reconstruction procedure, as 
well as the energy spectra, are discussed in ref. [18]. It is 
also of note that the data shown in panels (a) and (d) of 
figure]^ are the same as shown in figure 3 of ref. ini- 

It is clear from figure that the experimental data and 
the calculations agree in overall form, with the trend in 
the shape of the PADs with polarization well-reproduced 
by theory. This general behaviour is not surprising, since 
this sensitivity was inherent in the concept of obtaining 
the photoionization matrix elements via a fitting procedure 
from the 2D images recorded with different polarizations, 
but does indicate that the additional details observable in 
the tomographic data do not contradict the fit results, even 
at higher intensities which do affect the energy spectrum. 

The 2D images appear to show less satisfactory agree¬ 
ment but, since these (X, Y) plane projections include sum¬ 
mation over the Z-axis, this is perhaps unsurprising. In 
particular, the apparent increase in intensity of the band 
structures in panel (c), relative to the computational res¬ 
ults, is due to the additional (and incoherent) contribution 
from photoelectrons at different energies, due to the pro¬ 
jection of the broader spectrum onto the 2D plane, and 
which are not present in the computational results. The 
most significant differences are seen in panel (f), where the 
asymmetry in the (X, T) plane - the helicity of the distri¬ 
bution - is reduced relative to the computational results. 
This is likely due to a slight difference of the polarization 
ellipse relative to the calculations, as well as the summation 
over the broader spectrum (as mentioned above) which may 
wash-out fine details in the projection image. For the res¬ 
ults approaching circular polarization, panel (i), the agree¬ 
ment is better. In this case the contributions from the 
higher-energy AT feature are reduced (see ref. m), and 
the polarization state of the light may be slightly better 
matched to that assumed in the calculation. 

Overall these results indicate reasonably good agreement 
between the previously determined matrix elements and 
the tomographic data, but also indicate the problematic 
aspects of a qualitative comparison for these complex dis¬ 
tributions. In general such comparisons are worthwhile, 
but subject to perceptual bias which may be highly de¬ 
pendent on the type of data visualization used. Naturally 
a quantitative comparison is preferable, and is explored in 
the following section. 


B. Quantitative comparison 

To make a more careful comparison of the volumetric res¬ 
ults, were extracted from the data (see ref. [T8| for 
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Figure 7. Comparison of calculated & experimental tomographic results. Panels show (top) experimental and (middle) calculated 
3D distributions plotted as intensity iso-surfaces. The experimental data is plotted with a radial mask to show only the main part of 
the radial distribution. (Bottom) Experimental and calculated 2D projected images in the plane of polarization (X, Y). In this case 
no radial mask is included in the experimental results, and the images show an extended energy range relative to the computational 
results (the broader spectrum is due to the presence of Autler-Townes splitting in this case, see main text for details). Columns 
show the results for different laser polarization states, defined by a quarter-wave plate rotation angle (^A/ 4 ) in fhe experimental 
data, and spectral phase (j)y in the calculations. The states correspond to ellipticities e: ^0, 0.3 and 0.6 (see ref. [18] for further 
details). 


details). The /3 l,m(^) over the main feature can then be 
directly compared with the predicted (3^ ^ based on the 
fitted matrix elements. Since the fit results assume that 
the matrix elements are approximately constant over the 
feature, the experimental /5 l,m(^) were averaged over the 
FWHM of the main spectral feature to yield an energy- 
averaged value. These values are plotted in figure along 
with the calculation results. The range of the experimental 
/3L,M{k) are indicated by the error bars on the plot, indic¬ 
ating the spread of values over the spectrum. 

The agreement between the calculation and experimental 
results is generally very good, if not exact. The dominant 
terms, with L = 2, 4, 6 and M = 0, show excellent agree¬ 
ment and, aside from at = 0 , the experimental res¬ 
ults also show only a small spread of values. For the M 7 ^ 0 
terms, which generally have smaller magnitudes than the 
M = 0 terms except at large 0 ^, the agreement is generally 
good, but less so for the |M| = 2 terms. Here the trends 
with polarization are in agreement, but the exact values 
are shifted slightly from the calculations. As noted above, 
these small discrepancies may be due to slight differences 
in the laser polarization and frame rotations used in the 


calculations as compared with the experiment. 

Additionally, the experimental data indicated a signific¬ 
ant energy-dependence of the PADs away from the main 
spectral feature, small contributions from higher-order 
terms {L > 6 ) and symmetry breaking were also present. 
These effects are not accounted for by the net 3-photon 
model, and indicate the presence of additional complexit¬ 
ies to the light-matter interaction. These additional effects 
are, at this time, not well-understood beyond the clear re¬ 
quirement for ±m-state symmetry breaking, and for higher 
angular momentum states to be accessed. These observa¬ 
tions are discussed further in ref. m, and it is of note that 
the ability to resolve these additional effects via the quant¬ 
itative analysis of 3D photoelectron data is a significant 
outcome. 

Despite these additional, but intriguing, complexities, 
the major channels observed over the FWHM of the main 
spectral feature are seen to agree very well with the previ¬ 
ous analysis, based on 2D data recorded at lower intensities, 
overall providing a strong test of the accuracy of the ioniz¬ 
ation matrix elements determined in that case. The possib¬ 
ility of gaining a detailed understanding of the additional 
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Figure 8. Comparison of calculated ^lm (lines) and experi¬ 
mental results (crosses with error bars) as a function of polariz¬ 
ation state. Panels show L = 2, 4, 6, and line-styles and shade 
denote M. The error bars on the experimental results show 
the spread of ^L,M{k) obtained over the FWHM of the main 
spectral feature. 


effects observed in the 3D data, starting from the current 
3-photon model, and associated ionization matrix elements, 
remains an interesting proposition for future work. In the 
following section we explore some extensions to our treat¬ 
ment which may facilitate such understanding. 


V. ASSUMPTIONS, EXTENSIONS PHYSICAL 
CONSIDERATIONS 


In the above treatment, as applied in ref. m , some sim¬ 
plifications have been made for the specific case at hand, 
in order to facilitate the determination of the Rii(k) as de¬ 
tailed above (sect. 0- Here, in order to generalize this 
treatment further, we consider more carefully the assump¬ 
tions made and explore other extensions to the theory. 


Atoms 


The intra-pulse dynamics above implicitly assume that 
only the outermost electrons play a role in the intra-pulse 
dynamics, and that lower-lying bound states can be neg¬ 
lected. Furthermore, it assumes that the Ap manifold is the 
only unpopulated state which plays a role at the I-photon 
level. This is expected in this case because the As -^Ap 
transition carries significant oscillator strength, and is near 
resonant with the laser pulse. However, in general it is pos¬ 
sible that other states will play a role, particularly as the 
laser is tuned further from the As -^Ap line. This would 
result in a more complex TDSE (eqn. |^, with additional 
states appearing, and also necessitate a more careful treat¬ 
ment of the transition dipoles Pl/r to allow for variation 
in the transition amplitudes to different |n, /)-manifolds. 
The treatment of the 2-photon ionization would, simil¬ 
arly, increase in complexity with the addition of further 
initial/source |n,/)-manifolds, but would otherwise remain 
identical. 

In the case of atomic ionization with a structureless con¬ 
tinuum, the photoelectron energy spectrum can be treated 
somewhat directly as determined from the power spectral 
density of the laser pulse PU1I5U1I57] . Such treatment effect¬ 
ively introduces an additional time-dependent phase into 
eqn. In ref. m this phase is defined as with 

6uJe = ooe ^ oojp — uUp — 2 ujq , where the angular frequencies 
are related to the electron energy (huJe)^ Ihe ionization po¬ 
tential {hwip)^ the ionizing Ap state {hwp) and the photon 
energy {hw^). This phase will thus oscillate rapidly at the 
resultant difference frequency of these terms (effectively the 
difference between the total final state energy and the in¬ 
cident /input energy), and directly gives rise to a photoelec¬ 
tron energy spectrum dependent on the pulse properties, 
including its temporal duration and structure [38]. This 
dependence can be considered interferometrically, in the 
sense that the resultant (time-integrated) energy spectrum 
is the coherent temporal sum, hence contains interferences 
between all instantaneous momentum distributions; this is 
exactly analogous to the PADs considered as the coherent 
temporal sum of the instantaneous angular distributions 
(at a given energy). In the case of pulses intense enough 
to create significant Antler-Townes splitting in the photo¬ 
electron energy spectrum, this treatment could allow for 
a description of the changes in the PADs and symmetry 


breaking, as discussed in sect. |IVB| and correlated with 
the Antler-Townes doublet in the spectrum. This consid¬ 
eration is discussed further in ref. |I8| . 

In the most general case, where multiple, non-degenerate 
ionization pathways may be present, interferences may arise 
between ionizing transitions with very different angular 
structures. In the energy domain this effect has been in¬ 
vestigated by Elliott and co-workers in experiments utiliz¬ 
ing fundamental and second-harmonic light to create final- 
state interferences between different intermediates [39H4T] . 
Control over the relative phase of the two colours allowed 
for control over the resultant interferences uni . A sim¬ 
ilar concept was also employed to measure the phase of 
a bound-state [42]. Practically, this most general effect 
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could be included in our formalism by the inclusion of sets 
of ionization matrix elements correlated with the distinct 
sets of ionization pathways, where each set has a charac¬ 
teristic partial wave distribution (amplitudes and phases) 
and energy-dependent phase factor, and would result in the 
inclusion of interferences dependent on both geometric and 
energetic phase factors. Conceptually, this effect is inher¬ 
ent in the PADs arising from polarization shaped pulses, 
where the different ionization pathways correspond to dif¬ 
ferent intermediate angular momentum states, but in this 
case all levels are degenerate and the relevant phase shifts 
are purely geometric. 


B. Molecules 

In the case of molecular ionization the situation is more 
complex. In this case, the partitioning of the incident 
photon angular momentum to molecular rotations, as well 
as the outgoing photoelectron, requires a more involved 
treatment of the geometric terms, even in a single act¬ 
ive electron picture. Furthermore, one might expect that 
the continuum also contains structure due to population of 
different vibrational modes of the ion, although it is also 
possible that these states have little effect on the R{k) in¬ 
tegrals over a small energy range. This assumption form¬ 
ally means a separation of the ionization matrix element 
into electronic, vibrational (Franck-Condon) and rotational 
terms is possible, and that these terms are thus uncoupled. 
Effectively the electronic terms define the R{k)^ and the 
Franck-Condon factors an overall transition intensity en¬ 
velope - but one that does not affect the partial wave char¬ 
acter of the continuum. In some cases this approximation 
has been tested, and found to hold, but in other cases - par¬ 
ticularly when considering highly excited vibrational modes 
- one might expect this assumption to fail [43l |44] . 

In previous work we have investigated molecular ioniz¬ 
ation via vibrational and rotational state-resolved energy- 
domain experiments da ng na ES]. This work demon¬ 
strates the feasibility of performing such experiments, 
and illustrates the types of angular momentum coupling 
schemes required. Although the current work, incorporat¬ 
ing intra-pulse dynamics and polarization shaped pulses, 
has not yet been extended to molecular cases, such an ex¬ 
tension seems feasible based on these earlier studies, at least 
from the perspective of treating the ionization matrix ele¬ 
ments, and including the larger number of continuum l- 
waves required for molecular scattering problems. 

The intra-pulse dynamics in the molecular case may, 
however, be significantly more challenging. Clearly there 
are many more degrees of freedom to account for, and the 
potential for both nuclear and electronic wavepacket mo¬ 
tion during the pulse, as well as the possibility of dumping 
a lot of angular momentum into forming a rotational wave- 
packet during a strongly-coupled initial step (although the 
evolution of this wavepacket will ultimately be on a much 
slower timescale). The incorporation of such coupled rovi- 
bronic dynamics is, in practice, quite difficult due to the 
high-dimensionality of the problem. It is certainly not suffi¬ 


cient to perform a simple TDSE of coupled electronic states 
as employed herein, although the coupling of more com¬ 
plex wavepacket calculations with the ionization treatment 
herein would be feasible. A conceptually similar, although 
fully ab initio^ coupling of complex vibronic wavepackets 
with a full photoionization calculation has recently been 
presented for the triatomic molecule CS 2 EH; prior to this 
ab initio treatment a simpler dynamical model was com¬ 
bined with the relevant angular momentum coupling and 
gave good agreement with experimental results, although 
the treatment was only semi-quantitative and stopped short 
of extraction of the ionization matrix elements [48] . Eor di¬ 
atomics it is probable that a conceptual middle-ground is 
found, in which the required low-dimensionality wavepacket 
can be modelled via a simple TDSE treatment with enough 
accuracy to be of use. Eor polyatomics, the complexity 
of the wavepacket will be the deciding factor, depending 
directly on the number and type of states and couplings 
involved in a given case. In the most complex cases the 
problem may be best treated by fully ab initio calculations 
including photoionization, the results of which can be com¬ 
pared directly with experimental data at a high level, but 
in simpler wavepackets (few level and/or weakly coupled) 
a basic TDSE approach may be of sufficient accuracy to be 
useful for complete experiments. 

In sum, based on experience of similar problems in mo¬ 
lecular ionization in both the energy and time domains, it 
seems feasible that this time-domain multiplexing concept, 
employing multi-photon ionization schemes, and with the 
inclusion of intra-pulse dynamics, can also be applied suc¬ 
cessfully to (at least some) molecular photoionization prob¬ 
lems. 


C. Other regimes 

Other regimes are also of general interest in photoion¬ 
ization studies [T], in particular the strong-field (non- 
perturbative) regime. Very generally, the treatment presen¬ 
ted herein could be extended to this regime, and the issues 
associated closely follow the discussion above. In the non- 
perturbative case the intra-pulse dynamics become more 
complicated, since a single-active electron picture is no 
longer likely to be valid, and a static picture of the bound 
state energy-level structure also breaks down. Similarly, 
the scattering dynamics of the outgoing electron will also 
be time-dependent, since the scattering must now incor¬ 
porate the laser-induced part of the potential, not just the 
(static) atomic or molecular potential. 

In theory it is feasible to allow for these effects into 
the treatment presented herein, since it is already time- 
dependent and, as discussed for the atomic and molecular 
cases above, additional dynamical effects could be readily 
incorporated providing the numerics are tractable and ac¬ 
curate. However, the main issue in terms of determining 
the ionization matrix elements would be the large size of 
the set of matrix elements to be determined in a fully time- 
dependent treatment and the concomitant complexity of 
the fitting procedure if a set of ionization matrix elements 
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were required for each time-step. In such cases it may be 
possible to posit an effective functional dependence of the 
ionization dynamics on the laser field to mitigate this some¬ 
what, but one would have to more carefully consider exactly 
what kind of measurement would allow for a unique set of 
(fitted) matrix elements to be extracted from the (neces¬ 
sarily) time-integrated photoelectron image. 

Another regime is that of high-order light-matter coup¬ 
lings beyond the dipole approximation. In this case the 
ionization matrix elements contain higher-order angular 
momentum couplings, hence a more complex angular- 
momentum coupling scheme and a larger set of matrix 
elements must be determined, similar to the considera¬ 
tions for the molecular case. As for that case, there is no 
fundamental reason why such cases could not be treated 
within the theoretical framework presented here, although 
the feasibility of the fitting procedure would have to be as¬ 
sessed for any given case based on the size of the problem. 

Another interesting extension is to ionization time 
delays, since the the Wigner delay time is given by the 
energy-derivative of the scattering phase [49l [50]. Meas¬ 
urements of this phase, based on the concept of interfer¬ 
ing photoelectron wavepackets created with different en¬ 
ergies, have recently been demonstrated ED ED and are 
very similar to to concepts herein. In such measurements, 
above-threshold ionization creates electron wavepackets at 
different energies (i.e. multiple spectral features in the pho¬ 
toelectron spectrum), and further photo-absorption from a 
probe laser field can be used to interfere neighbouring wave- 
packets. This procedure results in side-band generation in 
the photoelectron spectrum, and the phase of the oscillation 
of these side-bands with respect to the probe field timing 
provides information on the relative phase of the photo¬ 
electron wavepackets - this is know as the RABBITT tech¬ 
nique [521 US] (this concept is somewhat analogous to the 
two-pulse photoelectron interferometry of ref. [30]). The 
difference between this concept and traditional “complete” 
photoionization experiments is that the total photoioniz¬ 
ation phase is measured over a broad energy spectrum in 
RABBITT measurements, as opposed to the measurement 
of the phases of the partial-waves at a single energy as 
discussed herein. By extending our technique to a broad 
energy range, e.g. via observation of multiple ATI features 
or the use of a broader bandwidth probe pulse, we would 
be able to obtain the partial-wave phases as a function of 
energy and, thus, determine the Wigner delay. Further¬ 
more, by obtaining the phases for all partial-waves, the 
angle-dependence of the Wigner delay in molecular ioniza¬ 
tion could also be investigated |54| . 


D. Maximum information measurements 
multiplexing 

In all cases discussed above, the main consideration is 
the feasibility of performing complete experiments for more 
complex, dynamical ionization schemes. Such applications 
will, naturally, be challenging, and require both a detailed 
theoretical understanding of the dynamics at hand and 


high-information experimental measurements. The major¬ 
ity of this work has focussed on assessing the results ob¬ 
tained for pure polarization states, in which there is no ad¬ 
ditional information gained from the coherent time-domain 
integration over the laser pulse, but for more complex 
cases the additional information content of polarization- 
multiplex measurements may be vital. Specifically, mul¬ 
tiplexing provides additional time-domain interferences in 
the PADs (see eqn. 12), with the result that time- 


integrated polarization-multiplexed measurements contain 
the information of multiple pure-state measurements. 

One particularly powerful aspect of using shaped pulses is 
the possibility of tailor-made pulses for metrology, designed 
to create or amplify specific interfering channels of interest. 
Conceptually this is identical to the use of shaped pulses for 
control EOlETl, however the pulses would be designed for 
the purposes of obtaining detailed information on specific 
ionization channels, rather than for the purposes of cre¬ 
ating a specific photoelectron distribution. The examples 
shown in figure and indicate how this concept oper¬ 
ates: there are different continuum populations created in 
the two cases, the time-domain structure is more complex 
in the shaped pulse case and, most generally, the pulse 
shape can be chosen to select certain ionization channels 
(within the constraints imposed by the dynamics of the ion¬ 
izing system). As well as polarization-shaping, the coherent 
time-domain treatment may also provide a way to probe 
additional intereferences due to effects such as intensity- 
dependent ionization phases. The presence of such effects 
at higher intensities has been determined from the beha¬ 
viour of the PADs over the Autler-Townes structure of th e 
photoelectron spectrum, as mentioned above (sect. |iVB| , 
but remains to be understood in detail. 

In all these cases, the PADs will usually be non- 
cyclindrically symmetric, so use of “maximum informa¬ 
tion measurements” utilizing 3D measurements and de¬ 
tailed analysis will also be required. The power of this 
approach has been touched on here, and is further explored 
in ref. EH]; recent work has also considered 3D photoelec¬ 
tron measurements in the context of photoelectron circular 
dichroism m 


VI. SUMMARY &: CONCLUSIONS 

In this work the validity of a fitting approach to complete 
photoionization experiments in the multiphoton regime, in¬ 
corporating intra-pulse dynamics, as initially reported in 
ref. [T7|, has been explored. The details of the fitting 
procedure, based on statistical sampling of the hyper¬ 
space and further testing and validation of the results, were 
outlined as a pragmatic fitting methodology. The results 
presented in ref. EH were discussed in detail, and compared 
both qualitatively and quantitatively with full 3D experi¬ 
mental photoelectron distributions. Finally, extension of 
this treatment to more complex ionization processes was 
discussed in general terms. 

This analysis indicated the validity of the results already 
presented, as well as insight into the practicalities of a prag- 
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matic fitting approach. Although this approach has yet to 
be tested beyond the use of pure polarization states, the 
use of polarization shaped pulses clearly offers an enhanced 
photoelectron metrology, with the possibility of controlling 
the information content via the pulse shape, as discussed 
in sect. (see also ref. [17]). The use of full 3D experi¬ 
mental measurements is another powerful aid to maximum 
information metrology, as indicated herein by comparison 
of the computational results with tomographically recon¬ 


structed experimental distributions (see also ref. (Hj). In 
general, we anticipate that the combination of these tools 
represents a powerful methodology for complete photoion¬ 
ization studies, or other research making use of ionization 
measurements. 
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